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Abstract 



We report results of ab initio calculation of the spin-rotational Hamilto- 
nian parameters including P- and P, T-odd terms for the BaF molecule. The 
ground state wave function of BaF molecule is found with the help of the 
Relativistic Effective Core Potential method followed by the restoration of 
molecular four-component spinors in the core region of barium in the frame- 
work of a non-variational procedure. Core polarization effects are included 
with the help of the atomic Many Body Perturbation Theory for Barium atom. 
For the hyperfine constants the accuracy of this method is about 5-10%. 
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a. Introduction. It is well known that possible P- and P, T-odd effects are strongly 
enhanced in heavy diatomic radicals (see, for example, In the molecular experiment 

with the TIF molecule stringent limits on the Schiff moment of the Tl nucleus and 
on the tensor constant of the electron-nucleus P, T-odd interaction were obtained. In the 
experiments with the polar diatomics with the unpaired electron one can search for the 
P, T-odd effects caused by the permanent electric dipole moment (EDM) of the electron 
de and by the scalar electron- nuclear P, T-odd interaction The most stringent limit 
on the electron EDM was obtained in the experiment with atomic Thallium (for the 
review of the theoretical predictions for dg see 0]). Heavy polar diatomic molecules provide 
enhancement of the electron EDM, which is several orders of magnitude larger, than in Tl. 
An experimental search for the EDM of the electron is now underway on the YbF molecule 
P]. The P-odd effects associated with the anapole moment of the nucleus are also strongly 
enhanced in diatomic radicals PJIO[]. 



The first calculations of the P, T-odd interactions in molecules were carried out for TIF 
molecule with the use of a "relativistic matching" of nonrelativistic one-configurational wave 



function 11 



Then a semi-empirical scheme [|T3,0] and ab initio method based on the Relativistic 



Effective Core Potential (RECP) calculation of the molecular wave function were de- 
veloped. The first RECP-based calculations of the P, T-odd spin-rotational Hamiltonian 
parameters for PbF and HgF molecules were carried out in the framework of the one- 
configurational approximation with minimal atomic basis sets, i.e. the correlation structure 
was not taken into account. In calculation of YbF molecule [0, a flexible atomic basis set 
was used and the correlation effects were considered within the Restricted Active Space SCF 
(RASSCF) method |T|Q. 



It was concluded in [|T^ that in order to perform more accurate calculations of the 
hyperfine and the P, T-odd constants, the spin-correlation of the unpaired electron with the 
outermost core shells 5s and 5p of ytterbium should be taken into account. Such correlations 
can be hardly efficiently considered within MC SCF-like methods because of the necessity 
to correlate too many electrons. 

Here we suggest to use an effective operator (EO) technique to account for the most 
important types of the core-valence correlations. EOs for the valence electrons are formed 
with the help of the atomic many body perturbation theory. This method allows to include 
correlations not only with the outermost core shells, but with all core electrons, which 
appears to be quite important for the hyperfine and P, T-odd interactions. The EO technique 
was recently developed for atoms |18| and proved to be very efficient for the calculations 



of the hyperfine structure of the heavy atoms \T^\. This technique is naturally and easily 
combined with the RECP method for the molecular calculations. As a result, a significant 
improvement of the accuracy is achieved. 

Below we report the results of application of this method to calculation of the BaF 
molecule. 

b. Spin-rotational Hamiltonian. Molecular spin-rotational degrees of freedom are de- 
scribed by the following spin-rotational Hamiltonian (see 0): 

ifsr = BN^ + 7SN - DeiiE + SAI 

+ WAkAn X S ■ I + {Wsks + Wdd,)Sn. (1) 
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In this expression N is the rotational angular momentum, B is the rotational constant, S 
and I are the spins of the electron and the Ba nucleus, n is the unit vector directed along the 
molecular axis from Ba to F. The spin-doubling constant 7 characterizes the spin-rotational 
interaction. and E are the molecular dipole moment and the external electric field. 
The axial tensor A describes magnetic hyperfine structure. It can be determined by two 
parameters: A = {A\\ + 2A±)/3 and A^ = {A\\ — A±)/3. The last three terms in account 
for the P- and P, T-odd effects. First of them describes interaction of the electron spin with 



the anapole moment of the nucleus /ca [|10|- The second one corresponds to the scalar P,T- 
odd electron-nucleus interaction with the dimensionless constant k^. The third one describes 
interaction of the electron EDM with the molecular field. Constant Wd characterizes an 
effective electric field on the unpaired electron. 

It is important to note that all P- and P, T-odd constants Wi mostly depend on the 
electron spin-density in the vicinity of the heavy nucleus. The same, of course, can be 
said about hyperfine constants A and A^. So, the comparison of the theoretical results for 
the hyperfine constants with the experiment is a good test for the accuracy of the whole 
calculation. 

c. RECP calculation of electronic wave function. The scheme of the RECP calculation 
for BaF molecule is very similar to that for YbF described in ||T^ (see also [|I^]) and below 
we will focus only on specific features of the present calculations. 

The Generalized RECP (GRECP) O (with the inner core ls\ . .]4:S^4:p^4:d^^ shells which 
were not included explicitly in the RECP calculations) was selected from a few other RECP 
variants for calculations of BaF because our test electronic structure calculations showed 
that it combined high accuracy with quite small computational expenses (see table |I| and 
the spectroscopic data below). 

Numerical pseudospinors derived from the GRECP/SCF calculations of some electronic 
configurations for Ba, Ba^ and Ba^+ were approximated by generally contracted s,p, d and 
/ gaussian functions forming (10, 8, 6, 2) [6, 5, 4, 2] basis set for bariumQ. For fluorine we 
used basis sets (14,9,4) — >• [6,5,2] and [4,3,3] from the ANO-I Library |1^. These basis 
sets proved to be sufficiently flexible to reproduce electronic structure in valence region of 
BaF as compared to other basis sets involved in our test SCF and RASSCF calculations. 

The RASSCF calculations of the spectroscopic constants were performed with the spin- 
Averaged part of the GRECP (AREP) and contribution of relatively small spin-orbit interac- 
tion (i.e. Effective Spin-Orbit Potential or ESOP as a part of GRECP) was estimated in the 
framework of the perturbation theory. The results of our AREP/RASSCF calculations with 
79558 configurations for the equilibrium distance and vibration constant [Rg. = 2.25 A, 
uOf, = 433 cm^^) are in a good agreement with the experimental data [^ (-Re = 2.16 A, 



^See [15j for details. Gaussian expansions for these pseudospinors, GRECP components and 



MO LCAO coefficients from BaF calculations can be found on http://www.qchem.pnpi.spb.ru 



^We used point group with (ai, foi, ^2) irreducible representations; 17 electrons were 
distributed on active orbitals within RAS 1=(3, 1,1,0), RAS 2=(3, 1,1,0) and RAS 3=(5,3,3,1) 
subspaces. 
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uOe = 469 cm~^). For the dipole moment we have obtained De = 2.93 D. 

d. Restoration of four- component spinor for valence electron. In order to evaluate 
matrix elements of the operators singular near nucleus of barium, we have performed 
GRECP/SCF and GRECP/RASSCF calculations of BaF where the pseudospinors corre- 
sponding to 5si/2, 5pi/2 and 5p3/2 shells were "frozen" with the help of the level-shift tech- 



nique (that is also known as Huzinaga-type ECP, see p2| and references therein). It was nec- 
essary to do because polarization of these shells were taken into account by means of EO tech- 
nique (see below) . Spin-orbit interaction was neglected for the explicitly treated electrons be- 
cause of its smallness (see [|T3|). Thus, only core molecular pseudoorbitals occupying mainly 
atomic Is, 2s and 2p orbitals of fluorine and the valence pseudoorbital of unpaired elec- 
tron (that is mainly 6s, 6p-hybridized orbital of barium) were varied. RASSCF calculations 
with 5284 configurations were performed for 11 electrons distributed in RAS 1=(2, 0,0,0), 
RAS 2=(2, 1,1,0) and RAS 3=(6,4,4,2) subspaces. 

The molecular relativistic spinor for the unpaired electron was constructed from the 
molecular pseudoorbital ip^ 

^^ = T.cm+T.cm"''=' + ---, (2) 

i i 

SO that the atomic s- and p-pseudoorbitals of barium in (^ were replaced by the unsmoothed 
four-component DF spinors derived for the same atomic configurations which were used in 
generation of basis s,p-pseudoorbitals. The MO LCAO coefficients were preserved after the 
RECP calculations. As the spin-orbit interaction for the unpaired electron is small, the 
"spin-averaged" valence atomic p-pseudoorbital was replaced by the linear combination of 
the corresponding spinors with j = / ± 1/2 (see [0,0 for details). 

e. Effective operators for valence electrons It is well known that the accuracy of the 
hyperfine structure calculations for heavy atoms is not high if core polarization effects are 



not taken into account. In ||2^ it was suggested, that correlations which are not included 
in the active space, can be treated with the help of the EO. The latter is constructed by 
means of the atomic many body perturbation theory (for the application of the perturbation 
theory to the calculations of the P, T- violation in atoms see, for example, |2^). The main 



advantage of this method is that there is no need to extend the active space to include core 
electrons. 

In [p3 it was supposed that EO is constructed in the active space which includes only 



few interacting levels. On the contrary, in [[18| , p!9[ it is suggested to use single EO for the 



whole (infinite dimensional) valence space. Thus, all correlations between valence electrons 
are treated explicitly, while EO accounts only for the core excitations. In this case, EO is 
energy dependent, but this dependence is weak if the energy gap between the core and the 
valence space is not too small. This makes EO method much more fiexible and allows to use 
one EO for different quantum systems, provided that they have the same core. In particular, 
it is possible to form EO for the atom (or ion) and then use it in a molecular calculation. 

Generally speaking, EO for the hyperfine interaction (as well as for any other one-electron 
operator) is no longer one-electron operator, even in the lowest order of the perturbation 
theory. On the other hand, the one-electron part of EO includes two most important cor- 
relation corrections and in many cases appears to be a very good approximation. The first 



4 



correction corresponds to the Random Phase Approximation (RPA), and the second one 
corresponds to the substitution of the Dirac-Fock orbitals by the Brueckner orbitals. 

To illustrate how EO works for the atomic barium, let us look at the hyperfine constant 
of the ^Pi(6s6p)-level of ^^^Ba. The two-electron multiconfigurational Dirac-Fock calculation 



gives A = 804 MHz which should be compared to the experimental value 1151 MHz. 
The two-electron configuration interaction calculation with RPA and Brueckner corrections 
included gives A = 1180 MHz. 

In this work we calculated EOs for the magnetic hyperfine interaction, for the EDM of the 
electron and for the anapole moment. Both RPA equations and Brueckner equations were 
solved for a finite basis set in the V^'"^ approximation (which means that SCF corresponds 
to Ba++), and matrix elements of the EOs were calculated. The basis set included Dirac- 
Fock orbitals for Is . . . 6s, 6p shells. In addition 7 — 21s, 7 — 21p, 5 — 20d and 4 — 15/ orbitals 



were formed in analogy to the basis set N2 of ||T8[. Molecular orbitals were reexpanded in 
this basis set to find matrix elements of EOs for the molecular wave function. 

/. Results. Expressions for the electronic matrix elements which correspond to the 
parameters A, A^ and Wi of the operator ([1|) can be found in 0. All radial integrals and 
atomic four-component spinors were calculated for the finite nucleus in a model of uniformly 
charged ball. 

Results for the parameters of the spin-rotational Hamiltonian are given in table |T|. There 
are two measurements of the hyperfine constants for ^^"^BaF [p6t^ . First of them was made 



for a matrix-isolated molecule and second was performed in a molecular beam. Results 



of these measurements were used in the semiempirical calculations [ll2|,p| of P- and P, T- 
odd parameters of the spin-rotational Hamiltonian. These calculations were based on the 
similarity between electronic matrix elements for the hyperfine structure interaction and for 
the P- and P, T-odd interactions. All of these operators mainly depend on the electron spin 
density in the vicinity of the nucleus. As a result, in a one-electron approximation parameters 
Wi are proportional to \/ AA^ |T^. Electronic correlations can break this proportionality. 

In table H we give results of the SCF and RASSCF calculations for 11 electrons with the 
restoration procedure described above. It is seen that in these calculations parameters A 
and Ad are significantly smaller than in experiments p6|p7| . On the next stage we used EOs 
to account for the core polarization effects. That led to the 50% growth for the constant A, 
while constant A^ increased by 130%. Our final numbers for the hyperfine constants are very 
close to the experiment (the difference being less than 5%) but differ more significantly 



from p7 |. 

Our SCF and RASSCF results for all three constants Wi are much smaller than results of 
the semiempirical calculations [0,0] • When core polarization effects are taken into account 
with the help of corresponding EOs, our values for Wa and Wa dramatically increase (at 
present we do not have RPA for the constant Ws)- There is a good agreement between our 
final value for Wd and that from the semiempirical calculation, but for the constant Wa, our 
result is noticeably smaller. 

It can be explained by the fact that proportionality between Wd and y/AAd holds within 
10% accuracy, but for the constant W\ deviation from proportionality reaches 30%. Almost 
half of this deviation is caused by the finite nuclear size corrections to radial integrals. 
Electron correlation corrections for both constants are about 15%. 

Two conclusions can be made from the results of this work. First, as it was suggested 
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in |]T5|, core polarization effects play very important role in calculations of parameters of 
the spin-rotational Hamiltonian for heavy diatomic radicals. Second, results of the ab initio 
calculations with core polarization included, are close to the results of the semiempirical 
calculations, correlation corrections being about 15%. The fact that two very different 
methods give similar results confirms that it is possible to make reliable calculations for 
such molecules. 
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TABLES 



TABLE L Excitation energies for low-lying states of Ba averaged over nonrelativistic configu- 
rations (finite difference SCF calculations). 





DF 


GRECP 




Transition 


Transition 


Absolute 


Relative 




energy (au) 


error (au) 


error (%) 


6s^ — > 6s^6p^ 


0.04813 


-0.00003 


0.06 




0.03942 


0.00010 


0.24 




0.15732 


-0.00002 


0.01 




0.24473 


0.00002 


0.01 




0.18742 


0.00017 


0.09 



TABLE II. Parameters of the spin-rotational Hamiltonian for BaF. 





A 






Wa 


Ws 




{MHz) 


{MHz) 


(1025 Hz ) 

e cm ' 


{KHz) 


{Hz) 


Exper .-I / Semiemp .'^ 


2326 


25 


-0.41 


240 


-13 


Exper.-II/Semiemp.*' 


2418 


17 


-0.35 


210 


-11 


SCF 


1457 


11 


-0.230 


111 


-6.1 


RASSCF 


1466 


11 


-0.224 


107 


-5.9 


SCF/EO 


2212 


26 


-0.375 


181 




RASSCF/EO 


2224 


24 


-0.364 


175 





^Hyperfine structure constants measured for matrix-isolated molecule [26| and semiempirical cal- 
culation of constants Wi based on this experiment H. 



'^Hyperfine structure constants measured for free molecule [27| and semiempirical calculation based 
on this experiment pi. 
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